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Abstract. Spectroscopy and V,I CCD photometry of the dwarf irregular galaxy SBS 1129+576 are presented for 
the first time. The CCD images reveal a chain of compact H n regions within the elongated low-surface-brightness 
(LSB) component of the galaxy. Star formation takes place mainly in two high-surface-brightness H n regions. 
The mean (V — I) colour of the LSB component in the surface brightness interval fiv between 23 and 26 mag 
arcsec -2 is relatively blue ~0. 56+0. 03 mag, as compared to the (V — J) ~0.9 - 1.0 for the majority of known 
dwarf irregular and blue compact dwarf (BCD) galaxies. Spectroscopy shows that the galaxy is among the most 
metal-deficient galaxies with an oxygen abundance 12 + log (O/H) = 7.36 ± 0.10 in the brightest H n region and 
7.48 ± 0.12 in the second brightest H II region, or 1/36 and 1/28 of the solar value*, respectively. H/3 and Ha 
emission lines and H<5 and H7 absorption lines are detected in a large part of the LSB component. We use two 
extinction-insensitive methods based on the equivalent widths of (1) emission and (2) absorption Balmer lines to 
put constraints on the age of the stellar populations in the galaxy. In addition, we use two extinction-dependent 
methods based on (3) the spectral energy distribution (SED) and (4) the (V — I) colour. Several scenarios of star 
formation were explored using all 4 methods. The observed properties of the LSB component can be reproduced 
by a stellar population forming continuously since 10 Gyr ago, provided that the star formation rate has increased 
during the last 100 Myr by a factor of 6 to 50 and no extinction is present. However, the observational properties of 
the LSB component in SBS 1129+576 can be reproduced equally well by continuous star formation which started 
not earlier than 100 Myr ago and stopped at 5 Myr, if some extinction is assumed. Hence, the ground-based 
spectroscopic and photometric observations are not sufficient for distinguishing between a young and an old age 
for SBS 1129+576. 
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- galaxies: individual (SBS 1129+576) 



1. Introduction 

SBS 1129+576 (a(J2000.0) = ll h 32 m 02!5, <S(J2000.0) = 
+57°22'45'.'7, Bicay et al. 2000) was discovered in the 
course of the Second Byurakan Survey (SBS) (Markarian 
& Stepanian 1983; Lipovetsky et al. 1988) as a galaxy with 
strong emission lines, weak continuum and ultraviolet ex- 
cess seen in a chain of H 11 regions embedded within an ex- 
tended blue low-surface-brightness (LSB) component. Up 
to now SBS 1129+576 has not been studied in detail. The 
low metallicity and relatively blue colour of its LSB com- 
ponent (this paper) make it one of the rare young dwarf 
galaxy candidates (Izotov & Thuan 1999). In the present 
paper the physical conditions and chemical abundances of 
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the ionized gas of SBS 1129+576 are studied for the first 
time. In addition, spectroscopic and V, I photometric data 
are used to study the properties of the unresolved stellar 
population in its bright H 11 regions and LSB component. 
Recently Thuan et al. (1999) derived for SBS 1129+576 
a redshift z = 0.00522 from single-dish H 1 21 cm ob- 
servations. After correction of the radial velocity for the 
Virgocentric infall motion, they derive a distance of D = 
26.3 Mpc, which we adopt. At this distance 1 arcsec cor- 
responds to a linear scale of 127 pc. 

The structure of the paper is as follows. In Sect. |3 we 
describe the observations and data reduction. The photo- 
metric properties of SBS 1129+576 are described in Sect. 
3. In Sect.0]we derive the chemical abundances in the two 
brightest H 11 regions. The properties of the stellar LSB 
population and its possible age range are discussed in Sect. 
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Fig. 1. (a) / image of the field around SBS 1129+576. The straight line shows the position of the slit during spec- 
troscopic observations. A pair of red background galaxies intersecting the LSB component ~15" northwards of the 
brightest region a (see panel b) is labeled bg. The ellipse delineates the central part of a probable uncatalogued back- 
ground galaxy cluster. Rectangles mark the position of some other uncatalogued galaxies in the field. These objects 
span a (V — I) colour range between 1.1 and 1.3 mag. (b) V image of SBS 1129+576. The overlayed contours corre- 
spond to surface brightness between 22.75 and 26.25 mag arcsec -2 in steps of 0.5 mag. The inset to the lower-right, 
computed from the V image after subtraction of a two-dimensional model for the LSB component, illustrates the 
spatial distribution of compact star-forming regions. The insets to the upper left and upper right show, respectively, 
close-up views in the V and / of the region depicted by rectangular brackets. This region shows a smooth morphology 
in the V, while the / image reveals two red ((V — 7)~1.6) background galaxies labeled #1 and #2. 



[SJ Finally, Sect. [B] summarises the main conclusions of this 
work. 

2. Observations and data reduction 

2.1. Photometric observations and data reduction 

Direct images of SBS 1129+576 in V and I (Fig.[TJ) were 
acquired with the Kitt Peak 2.1m telescope 1 on April 
19, 1999, during a photometric night. The telescope was 
equipped with a Tektronix 1024x1024 CCD detector op- 
erating at a gain of 3e~ ADU -1 , giving an instrumental 
scale of 0'.'305 pixel -1 and field of view of 5'. The total 
exposures of 20 and 30 min in V and /, respectively, were 
split into four subexposures each being slightly offset with 
respect to each other for removal of cosmic particle hits 
and bad pixels. The point spread function in V and / 

1 Kitt Peak National Observatory (KPNO) is operated by 
the Association of Universities for Research in Astronomy 
(AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF). 



were respectively 1'.'06 and l'/17 FWHM. Bias and flat- 
field images were obtained at the beginning and end of 
night. Calibration was achieved by observing 4 standard 
fields from Landolt l|1992fl at 3-4 different airmasses dur- 
ing the night. Our calibration uncertainties are estimated 
to be 0.01-0.02 mag in each of the V and I bands. The 
data reduction, including bias subtraction, removal of cos- 
mic particle hits, flat-field correction and absolute flux 
calibration was made using IRAF 2 . 

2.2. Spectroscopic observations and data reduction 

The spectroscopic observations were carried out on 
19 June, 1999, with the Kitt Peak 4m Mayall tele- 
scope and Ritchey-Chretien spectrograph with the T2KB 
2048 x 2048 CCD detector, with the slit at P.A. = 171°, 
centered on the brightest star-forming region and ex- 

2 IRAF is the Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the AURA under cooperative agreement 
with the NSF. 
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tending along the elongated body of the galaxy (roughly 
aligned with the major axis; see Fig.[TJi). A 2" x 300" slit 
with the KPC-10A grating in first order and a GG 375 
order separation filter was used. The spatial scale along 
the slit was 0'.'69 pixel -1 and the spectral resolution ~7 A 
(FWHM). The spectra were obtained at an airmass 1.33 
and in a total exposure of 60 minutes, which was broken 
up into 3 subexposures. No correction for atmospheric re- 
fraction was made because the slit was oriented with a 
P.A. close to the parallactic angle. Two Kitt Peak spec- 
trophotometric standard stars were observed for flux cali- 
bration. For wavelength calibration, spectra of a He-Ne-Ar 
comparison lamp were taken after each exposure. 

The data reduction was made with the IRAF software 
package. This includes bias-subtraction, flat-field correc- 
tion, cosmic-ray removal, wavelength calibration, night 
sky subtraction, correction for atmospheric extinction and 
absolute flux calibration of the two-dimensional spectra. 

One-dimensional spectra for abundance determination 
in the two brightest H II regions a and b (Figs.^3 and 01 
were extracted from the two-dimensional spectrum within 
large apertures of 2" x 5". Some additional spectra of 
regions a and b within smaller apertures 2" x l'.'l, 2" x 
1"4 and 2" x 2'/l were also extracted. 

In addition we extracted spectra showing hydrogen 
B aimer absorption lines for five regions along the ma- 
jor axis of the galaxy to study the stellar population of 
the LSB component. The locations of selected regions, de- 
noted 1 to 5, relative to region a are given in Tables 0] and 
03 The spatial extent of these regions along the slit were 
6"2, 5"5, 4"8, 4"8 and 4"8, respectively. 

3. Photometric analysis 

3.1. Morphology, environment and colour distribution 

SBS 1129+576 appears as an oblate dwarf {My ~ -15.9 
mag) galaxy with a projected axial ratio of 4 at 25 V mag 
arcsec -2 (Fig. ^ HJ- Star formation activity takes place 
in the central parts of the galaxy within an elongated, 
moderately blue region with a size of ~ 2.5 kpc. The V 
magnitudes of regions a and b are respectively f» 19.0 mag 
and 20.0 mag, corresponding to absolute V magnitudes - 
13.1 mag and -12.1 mag. The (V — I) colours of ~ 0.1 
mag and ^0.3 mag, respectively, of regions a and b are 
relatively blue. The (V — I) colour of the LSB component 
shows a smooth change from ~ 0.44 mag in the immedi- 
ate vicinity of region a to an average colour of ~ 0.56 ± 
0.03 mag in the outer part of the galaxy (see Fig.^). The 
(V — I) colour map (Fig. 0) reveals for surface brightness 
levels fainter than 24 V mag arcsec -2 , a featureless and 
relatively constant colour over the whole LSB component 
except for a strikingly red ((V — I) ~ 1.6 mag) region lo- 
cated ~ 15" north of region a (region bg in Figs. and [2]). 
The local colour excess observed in region bg is due to two 
background galaxies seen in the I image only. A close-up 
view of this region in the V and I is shown in the upper 
insets of Fig. \iy>. A potential slight overestimate of the 
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Fig. 2. (V - I) map of SBS 1129+576 displayed in the 
colour range between +0.35 to +0.8 mag. The overlayed 
contours are from 21.1 to 25.0 V mag arcsec -2 in steps 
of 0.3 mag. The region, depicted by the ellipse, shows a 
strikingly redder (V — I) colour (~ 1.6 mag) than the 
average colour ~ 0.56 mag of the LSB component. This 
local maximum in the (V — 7) colour is to be attributed 
to background galaxies labeled bg in Fig. and in the 
upper-right inset of Fig. ^p. 



(V-I) colour of the LSB component in SBS 1129+576 as 
a result of the superposition of background sources at dif- 
ferent locations is likely given the numerous red ((V — I) 
> 1.2 mag) faint (mj 20.5 mag) sources in the field of 
the galaxy. Figure^ shows that SBS 1129+576 is located 
in front of a probable background cluster of galaxies cen- 
tered ~ 50" west of region a. This cluster with the central 
part delineated by the ellipse as well as several red galax- 
ies indicated by rectangles in Fig. are not catalogued 
in the NASA/IPAC Extragalactic Database (NED). 

3.2. Surface photometry 

Surface brightness profiles (SBPs) of SBS 1129+576 have 
been computed following the methods i through iii de- 
scribed in Papaderos et al. I|1996afl . Briefly, the photo- 
metric radius R* = (A(fi)/iT) - 5 corresponding to the sur- 
face brightness level fi is computed from the area A(p) 



4 N. G. Guseva et al. 



Table 1. Structural properties of the starburst and LSB components of SBS 1129+576 a . 



Band 


A*E,0 
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E-25 
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mtot 






mag arcsec 


pc 


kpc 


mag 


kpc 


mag 


mag 


mag 


mag 


kpc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


V 


21.01±0.02 


436±4 


0.73 


18.27 


1.59 


16.54 


16.43 


16.23±0.01 


16.22 


0.65,1-17 


I 


20.43±0.03 


433±4 


0.78 


18.25 


1.82 


15.92 


15.86 


15.73±0.02 


15.72 


0.69,1.20 



a The tabulated values have not been corrected for interstellar extinction or inclination. 
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Fig. 3. (a) Surface brightness profiles (SBP) of SBS 1129+576 in the V and / derived using the method iii. The 
modeled surface brightness distribution of the LSB component in the V assuming a modified exponential distribution 
with (6,g)=(1.2,0.4) is shown by the thick-grey curve. Open circles show the surface brightness distribution of the light 
in excess to the modeled LSB profile and is to be attributed to the compact star-forming regions along the major axis 
of the galaxy (Fig.GJs). The effective radius r e ff of the V SBP and the isophotal radii P25 and E25 of, respectively, the 
starburst and the LSB component at 25 V mag arcsec -2 are indicated. The dashed line shows a linear fit to the / band 
profile for R* >8", extrapolated to R* — 0. (b) Circles show the (V — I) colour profile computed by subtraction of the 
SBPs displayed in the left panel. Labels attached to the upper axis indicate the V surface brightness corresponding to 
the respective photometric radius R*. The (V — I) colour distribution (thick-grey line) is obtained by subtraction of 
the fits in V and / shown in (a) . It closely matches the observed colour distribution and implies a practically constant 
(V — I) colour of w 0.56 mag over a range of 4 magnitudes in surface brightness. 



of a galaxy in arcsec 2 , as derived through ellipse fitting or 
computation of a line-integral along an isophote (methods 
i and ii) or summation of all pixels inside a polygonal aper- 
ture with a surface brightness brighter than fj, (method 
iii). Essentially, these techniques trace the growth of the 
isophotal size of a galaxy with decreasing intensity /. They 
require no choice of a "geometrical center" of a galaxy 
and insure that the photometric radius R* is a monotonic 
function of ji. Evidently, in order to derive SBPs as de- 
scribed above, one has to keep track of the morphology of 
a BCD throughout its intensity span, i.e. in general to be 
able to interpolate an isophote down to the faintest mea- 
sured level /z of a SBP. This allows to visually check for 
and screen-out foreground and background sources in the 
periphery of the galaxy, thus to make sure that source con- 
fusion does not affect the SBP slope at faint intensity lev- 
els. This task is more difficult to achieve when computing 



SBPs based on photon statistics inside circular or ellipti- 
cal annuli, extending out to a user-defined maximal radius 
''max- Especially for BCDs, SBPs derived with the latter 
methods may considerably vary, depending on subjective 
assumptions on the r max or the "center" of a galaxy. 

As a check for consistency, we also computed SBPs 
using method iv in Papaderos et al. (2002). This technique 
is based on the calculation of photon statistics for a series 
of masks of arbitrary (generally irregular) shape, mapping 
equidistant logarithmic intensity intervals between I m i n 
and / m ax- As in methods i through iii, method iv does not 
require the choice of a "geometrical center" and accounts 
adequately for the large morphological variation of a BCD, 
with a typically smooth LSB part and an irregular star- 
forming component. 

We derived V and / SBPs from the total light of the 
galaxy, except for region bg (shown by the ellipse in Fig. 
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|2J) which has been replaced by a two-dimensional fit to the 
adjacent LSB emission. Both SBPs in Fig. [3]} are derived 
using the method iii. They show an exponential intensity 
decrease for radii R* <; 8" with a scale length a «3'.'4. 
That the exponential model yields a reasonable approxi- 
mation to the LSB emission is also indicated from fitting 
a Sersic 0968) profile 



I(R*) = I exp 



R* 
T 



(1) 



(see also Caon et al. 119931 Cellone et al. 119941 Papaderos 
et al. H996a|l to the radius interval R* >8". The exponent 
7/ obtained this way, respectively 1.2 and 1.1 in V and I, 
is marginally larger than the value rj — 1, corresponding 
to the exponential law. 

However, a closer inspection of the I band SBP (Fig. 
|3Ji) shows that an inwards extrapolation of the exponential 
law fitted to the LSB component (see dashed line, over- 
layed with the I profile) predicts for small radii (i?*^2") 
a slightly higher intensity than the one observed. Such a 
pure exponential LSB model would imply that the star- 
forming component in SBS 1129+576 provides no more 
than 5% of the total / emission and that its contribu- 
tion decreases rapidly at small radii. This suggests that 
the stellar LSB cmision of SBS 1129+576 is best approx- 
imated by an exponential profile with a flattening in its 
inner part. Note that such SBPs have been frequently ob- 
served in dwarf ellipticals (e.g., Binggeli & Cameron 1991 ), 
dwarf irregulars (Ronnback & Bergvall 1994; Patterson 
& ThuanUHMI Makarova van ZeeEUOl) and a few 

blue compact dwarf (BCD) galaxies (e.g., Papaderos et al. 
Il99tjal Vennik et al. [TUM I2UUU1 Telles & Terlevich ITM71 
Guseva et al. lMTTl Fricke et al. 12003) ■ SBP s of th is kind, 
classified "type V" in Binggeli & Cameron (|1991fl , can be 
approximated by, e.g., the modifed exponential distribu- 
tion proposed in Papaderos et al. 0996a): 

I (If) = I E ,oexp (—} {1 - qcxp(-P 3 (R*))} (2) 



with 



R* 1 — q 



(3) 



a q 

The empirical fitting function (Eq. |2J) flattens with re- 
spect to the exponential law inside of a cutoff radius ba, 
and attains at R*=0" an intensity given by the relative 
depression parameter q — A7/Jeo < L An advantage of 
Eq. [21 is that its exponential part, depending on Ie$ and 
a only, can be readily constrained from linear fits to the 
outer exponential part of a "type V" SBP. 

In order to disentangle the intensity distribution of the 
LSB component and better constrain the depression pa- 
rameters b and q in Eq. [5] we follow the approach adopted 
in Guseva et al. ( 2001). We first subtracted compact (di- 
ameter 554") high-surface-brightness regions in the inner 
part of the galaxy and then rederived SBPs from the resid- 
ual underlying LSB emission. Fitting Eq. [3 to the result- 
ing profiles yields V and I depression parameters (b,q) — 



(1.2,0.4) and an exponential scale length a ~ 430 pc. In 
Fig.|3Ji we show for the V SBP the modeled surface bright- 
ness distribution of the LSB component according to Eq. 
13 and the emission in excess of the model with the thick- 
grey curve and open circles, respectively. The dashed line 
shows a linear fit to the I band profile for R* >8", extrap- 
olated to R* =0. The excess emission is due to the chain 
of compact star-forming regions along the major axis of 
the galaxy. 

Table ^ summarizes the derived photometric quanti- 
ties. Cols. 2 and 3 give, respectively, the central surface 
brightness /ie,o and scale length a of the LSB component 
as obtained from linear fits to the SBPs for R* > 8" and 
weighted by the photometric uncertainty of each point. 
These quantities correspond to the values one would ob- 
tain from extrapolation of the exponential slope observed 
in the outlying regions of the galaxy to R* — 0. Cols. 4 
through 8 list quantities obtained from profile decompo- 
sition where the intensity distribution of the LSB compo- 
nent was modeled by the modified exponential distribu- 
tion (Eq. 0). Cols. 4 and 6 give the radial extent P25 and 
E25 of the starburst and LSB components, respectively, 
both determined at a surface brightness 25 mag arcsec~ 2 . 
The apparent magnitudes of both components within P25 
and E25 are listed in cols. 5 and 7, respectively. Col. 8 
gives the apparent magnitude of the LSB component in 
each band within a photometric radius of 18", as derived 
from integration of the modeled distribution (Eq. |2J) . The 
total magnitudes of the galaxy, as inferred from SBP inte- 
gration out to the same radius and by integrating the flux 
within a polygonal aperture are listed in cols. 9 and 10, 
respectively. Col. 11 gives the effective radius r c g and the 
radius r§o which encircles 80% of the galaxy's total flux. 

From Tabled it is evident that the integrated emission 
of the starburst component including the two brightest 
regions a and b contributes only ~ 17% of the V light 
of SBS 1129+576 within its 25 V mag arcsec~ 2 isophote. 
This is a factor of 3 lower than the average of 50% derived 
for BCDs in the B band by Papaderos et al. I)1996b|) and 
Salzer & Norton <|T39"9|l . 

The (V - I) colour profile of SBS 1129+576, derived 
from subtraction of the corresponding SBPs, is shown 
in Fig. |5Jd. Its behaviour is similar to that in BCDs 
(Papaderos et al. I1996a|l and compact irregular dwarf 
galaxies (Patterson & Thuan 115561 van Zee |2"UUU|) with 
a blue colour in the inner part of the galaxy and a red- 
der, relatively constant colour at larger radii. The colour in 
SBS 1129+576 increases gradually from (V — I) 0.2 mag 
for radii R* ^ 2" to ~ 0.5 mag at R* ^ 5" and remains 
practically constant at ~ 0.56 ± 0.03 mag in the outer 
part of the galaxy. The observed colour is in accord with 
the one resulting from subtraction of the apparent magni- 
tudes of the modeled distributions for the LSB component 
(Tabled col. 8) which has an average value of 0.57 mag. 
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Fig. 4. The KPNO 4m spectra of regions a and b with the emission lines identified. The lower spectra in (a) and (b) 
are the observed spectra downscaled by a factor of 15. 



4. Chemical abundances 

In this Section we analyze the element abundances in SBS 
1129+576 based on spectroscopic observations of the two 
brightest H n regions a and b. The spectra of these star- 
forming regions (Fig.^l are characterised by strong nebu- 
lar emission lines superposed on stellar Balmer absorption 
lines. The latter are also seen along the slit in the LSB 
component. 

The emission line fluxes were measured using Gaussian 
profile fitting. The errors of the line flux measurements in- 
clude the errors in the fitting of profiles and those in the 
placement of the continuum. They have been propagated 
in the calculations of the elemental abundance errors. The 
observed (F(X)) and corrected (/(A)) emission line fluxes 



relative to the H/3 fluxes, the equivalent widths EW of the 
emission lines, the observed fluxes of H/3, and the equiva- 
lent widths of the hydrogen absorption lines are listed in 
Table □ 

The Ha-to-H/3 flux ratios in both H n regions are lower 
than the theoretical value (e.g., Brocklehurst 1971). This 
is likely not due to data reduction problems as the Ha-to- 
H/3 flux ratios in other galaxies observed during the same 
night are greater than the theoretical ones. Therefore, 
an extinction coefficient C(H/3)=0 was assumed for these 
H II regions and the emission-line fluxes were corrected for 
Balmer line absorption only. 

Some diagnostic lines were studied to check for possi- 
ble deviations of the H II region emission in SBS 1129+576 
from the predictions of photoionization models. For this 
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Table 2. Observed (F(X)) and corrected (/(A)) fluxes and equivalent widths (EW) of emission lines in regions a and 
b. 



Ao(A) Ion 


region a 




region b 




F(A)/F(H/3) /(A)//(H/3) 


£W(A) 


F(A)/F(H/3) I(X)/I(Hf3) 


EW(A) 



3727 


[O ii] 


1. 


.408 


±0. 


.035 


1 


.346 


±0 


.037 


51. 


.0 


±0. 


.8 


1 


.875 


±0. 


.061 


1 


.758 


±0. 


.065 


38. 


.8 


±0 


.7 


3835 


H9 





.030 


±0. 


.013 





.084 


±0 


.048 


1. 


.0 


±0. 


.4 


























3868 


[Nc in] 


0. 


.169 


±0. 


.015 





.161 


±0 


.015 


4 


.5 


±0. 


.4 





.313 


±0. 


.026 





.293 


±0. 


.026 


5. 


.3 


±0 


.4 


3889 


H8 + He i 





.128 


±0. 


.014 





.179 


±0 


.026 


4 


.2 


±0 


.5 





.147 


±0. 


.025 





.240 


±0 


.051 


3 


.1 


±0 


.5 


3968 


[Ne in] +H7 


0. 


.122 


±0. 


.012 





.177 


±0 


.024 


3. 


.7 


±0. 


.4 





.158 


±0. 


.027 





.249 


±0. 


.052 


3. 


.3 


±0 


.6 


4101 


H<5 





.183 


±0. 


.013 





.235 


±0 


.023 


5. 


.5 


±0 


.4 





.180 


±0. 


.020 





.252 


±0. 


.044 


3 


.5 


±0 


.4 


4340 


H 7 





.457 


±0. 


.018 





.490 


±0 


.024 


15. 


.6 


±0. 


.5 





.434 


±0. 


.026 


0. 


.483 


±0. 


.038 


9. 


.1 


±0 


.5 


4363 


[O in] 


0. 


.061 


±0. 


.014 





.058 


±0 


.014 


1 


.7 


±0. 


.4 





.117 


±0. 


.023 





.110 


±0. 


.023 


2 


.2 


±0 


.4 


4471 


He I 





.033 


±0. 


.014 





.032 


±0 


.014 


1. 


.0 


±0. 


.4 





.042 


±0. 


.019 





.038 


±0. 


.019 





.8 


±0 


.3 


4861 


H/3 


1. 


.000 


±0. 


.027 


1 


.000 


±0 


.030 


41 


.0 


±0 


.8 


1 


.000 


±0. 


.040 


1 


.000 


±0 


.045 


25. 


.4 


±0 


.8 


4959 


[O in] 


0. 


.666 


±0. 


.022 





.637 


±0 


.022 


24 


.4 


±0. 


.6 


1 


.173 


±0. 


.044 


1 


.099 


±0. 


.044 


26. 


.6 


±0 


.7 


5007 


[O in] 


1. 


.975 


±0. 


.045 


1 


.888 


±0 


.045 


73 


.8 


±0. 


.9 


3 


.252 


±0. 


.099 


3. 


.049 


±0. 


.099 


76. 


.9 


±1 


.0 


5876 


He i 


0. 


.086 


±0. 


.010 





.082 


±0 


.010 


4. 


.5 


±0. 


.5 





.093 


±0. 


.017 


0. 


.087 


±0. 


.017 


3. 


.0 


±0 


.6 


6300 


[O i] 


0. 


.048 


±0. 


.010 





.046 


±0 


.010 


3 


.0 


±0 


.6 





.075 


±0. 


.017 





.069 


±0 


.017 


3 


.1 


±0 


.7 


6563 


Ha 


2 


.692 


±0. 


.058 


2 


.599 


±0 


.064 


186. 


.3 


±1 


.8 


2 


.662 


±0. 


.079 


2 


.526 


±0 


.088 


137. 


3 


±1 


.6 


6584 


[Nil] 


0. 


.041 


±0. 


.008 





.039 


±0 


.008 


2. 


.8 


±0 


.5 





.071 


±0. 


.013 





.066 


±0 


.013 


3 


.1 


±0 


.6 


6678 


He i 


0. 


.030 


±0. 


.008 





.029 


±0 


.008 


2. 


.1 


±0. 


.6 





.026 


±0. 


.015 


0. 


.024 


±0. 


.015 


1 


.1 


±0 


.7 


6717 


[S ii] 


0. 


.142 


±0. 


.010 





.136 


±0 


.010 


10. 


.1 


±0. 


.7 





.168 


±0. 


.018 


0. 


.155 


±0. 


.018 


7. 


.5 


±0 


.8 


6731 


[S ii] 


0. 


.103 


±0. 


.010 





.098 


±0 


.010 


7. 


.3 


±0. 


.7 





.114 


±0. 


.015 


0. 


.106 


±0 


.015 


5. 


.1 


±0 


.6 



C(R(3) dex 
F(H/3) a 
EW(abs) A 



0.000±0.028 
0.50±0.01 
1.9±0.4 



0.000±0.041 
0.23±0.01 
1.7±0.4 



Mn units 10 14 erg s 'cm 2 . 



purpose, data for the emission lines were collected for 
the galaxies from Izotov et al. (|19941 11997(1 . Thuan ct 
al. ifT^B) . Iz otov fc Thuan (ESHJ), Frickeet al. <PUT|1 . 
Noeske et al. (2000), and Guseva et al. (2001). In total we 
use data for 46 star- forming galaxies: 11 lowest-metallicity 
galaxies with an oxygen abundance 12 + log(0/H) < 7.6 
and 35 higher-mctallicity galaxies with an oxygen abun- 
dance 12 + log(0/H) > 7.9. Various emission line flux 
ratios are shown in Fig. [S] 

In this figure, galaxies with low and high oxygen abun- 
dances split into two sequences, shown by filled and open 
circles, respectively. This separation is in overall agree- 
ment with photoionization H II region models. However, 
regions a and b of SBS 1129+576 (filled triangles in Fig.|3J) 
having lower ionization parameters, lie outside the main 
location of low-metallicity galaxies. The line ratios [S 
ll]A6717+6731/Ha and [N n]A6584/Ha in SBS 1129+576 
are ~ 3 times and [O i]A6300/Ha ~ 5 times higher than 
the ones in the sample of the 11 lowest-metallicity galaxies. 
The locations of regions a and b in Fig. [5] also significantly 
deviate from the model predictions calculated with the 
CLOUDY code (Ferland Ferland et aL Eggj) for a 

heavy element mass fraction Z = Zq /36 derived for region 
a and for ionizing stars with effective temperatures of T e ff 
= 45 000 K (thin solid lines) and T eS = 50 000 K (thick 
solid lines). By dot-dashed and dashed lines we also show 
the model predictions for T c g = 50 000 K and two heavy 
element mass fractions of Z@/20 and Zq/5, respectively. 



The discrepancies between photoionization models and 
the observed line intensity ratios (essentially the high [O 

I] A6300/H/5 and [S n]A6717+6731/H/3) are usually ex- 
plained by some contribution of shock waves. However, 
Stasihska & Izotov ( 2003 ) from the analysis of a large sam- 
ple (~ 400) of H II galaxies proposed another model to ex- 
plain the enhancement of these lines in low-metallicity H 
II regions by invoking chemical inhomogeneities and self- 
enrichment by the heavy elements. Their models take into 
account the time evolution of an ionizing cluster in the 
simple case of an expanding bubble. The models predict 
an increase of the [O i]A6300/H/3 ratio in better agree- 
ment with observations of low-metallicity H II regions. 
Nevertheless, the hypothesis of shock heating may also be 
considered as an alternative explanation. 

A two-zone photoionized H II region model has been 
assumed for the abundance determination. The elec- 
tron temperature T (0 in) in the high-ionization re- 
gion has been derived from the observed flux ratio [O 
m]A4363/(A4959+5007), using a five-level atom model 
fAUer lrggl)! with atomic data from Mendoza ljlUg5|l . The 
electron temperature T c (0 ii) in the low- ionization region 
has been obtained using the empirical relation between 
T (0 ii) and T e (0 in) from the H II region photoioniza- 
tion models by Stasihska flMty . The [S n]A6717/A6731 
ratio was used to derive the electron number density N C (S 

II) . 

The abundances of +2 , Ne +2 and He + were derived 
applying the electron temperature T c (Oin). The electron 
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-1.5 -1 -0.5 0.5 



log([OII]/[OIII]) 

Fig. 5. Dependence on a logarithmic scale of the line flux 
ratios (a) log([N n]A6584/Ha) and (b) log([0 i]A6300/Ha) 
on the log([0 n]A3727/[0 m]A5007) ratio. Data for galax- 
ies with 12 + log(0/H) < 7.6 and 12+log(0/H) > 7.9 are 
shown by filled and open circles, respectively. The posi- 
tions of the brightest regions a and b in SBS 1129+576 
are shown by triangles. Thin and thick solid lines rep- 
resent theoretical predictions from photoionization mod- 
els calculated with the CLOUDY code for ionizing stars 
with effective temperatures T e g = 45 000 K and T c g = 
50 000 K and a heavy element mass fraction Z = Zq /36. 
Dot-dashed and dashed lines correspond to models with 
T e s = 50 000 K and a heavy element mass fractions Z = 
Zq/20 and Zq/5, respectively. The ionization parameter 
decreases from log U — -2.0 (left) to log U = -3.0 (right). 

temperature T c (On) is adopted for the + and N + ionic 
abundance determination. 

The total oxygen abundance is the sum of the + and 
+2 abundances. The total abundances of other heavy ele- 
ments were derived using ionization correction factors fol- 
lowing Izotov et al. lfT^llrW7jl and Thuan et al. ijlMS) . 
The ionic and heavy element abundances for regions a and 
b together with electron temperatures and electron num- 
ber densities are given in Table [3] along with the adopted 
ionization correction factors (ICF). 

The oxygen abundances 12 + log(0/H) = 7.36 ± 0.10 
(Z=Z /36) and 7.48 ± 0.12 (Z=Z /28) for regions a and 
b, respectively, are slightly different but still consistent 



within the la errors. Note that the oxygen abundance in 
region a may be underestimated due to the additional con- 
tribution of shock enhancement of the [O iii]A4363 emis- 
sion line. Such an effect can be present in region 6, but it 
is expected to be much larger in region a with the much 
weaker [O iii]A4363 emission line (Table 

The abundance ratios of nitrogen and neon to oxygen 
(Table |3J are in agreement with those derived in other 
lowest-metallicity dwarf galaxies (Izotov & Thuan 1999). 

The low metallicity of SBS 1129+576 makes it poten- 
tially suitable for He abundance determination. For this, 
we use the three strongest He I A4471, A5876 and A6678 
emission lines. The helium abundance derived from the 
corrected fluxes of all observed He I emission lines is shown 
in Table |3 The mean values for the He mass fraction of 
Y = 0.220 ± 0.027 (region a) and Y = 0.202 ± 0.042 (re- 
gion b), though formally consistent, are significantly lower 
than the primordial He mass fraction Y p = 0.244 ± 0.002, 
derived by extrapolating the Y vs O/H linear regression 
to O/H = for the sample of 45 low-metallicity dwarf 
galaxies (Izotov & Thuan lT^Sj l. or to Y p = 0.245 ± 0.002 
derived for the most metal-deficient BCDs I Zw 18 and 
SBS 0335-052 (Izotov et al. lTM7|> . 

The likely source of lower values of Y in SBS 1129+576 
is underlying stellar He I absorption superposed on the 
He I emission lines. Comparing the equivalent widths of 
He I emission lines (Tabled} with those predicted by pop- 
ulation synthesis models for He I absorption lines (e.g., 
Gonzalez Delgado et al. 1999) we conclude that underly- 
ing stellar absorption can decrease the fluxes of the He I 
emission lines by as much as ~ 10 - 20% for He I A5876 
(assuming the equivalent widths of He I A5876 and He I 
A4471 absorption lines to be similar) and more for other 
He I lines. Hence, despite its low metallicity SBS 1129+576 
is not a good candidate for primordial helium abundance 
determination because of the large effect of the underlying 
stellar absorption. 

5. Age of the underlying stellar population 

The low metallicity and relatively blue (V — I) colour of 
SBS 1129+576 make this object a good young galaxy can- 
didate. In this section we consider the properties of the 
stellar populations in the galaxy and discuss its evolution- 
ary status. The observed properties of stellar populations 
are dependent on the metallicity and star formation his- 
tory. They also can be influenced by interstellar extinction 
and emission of the ionized gas. Therefore, we use all avail- 
able spectroscopic and photometric observational data on 
5 LSB regions to put consistent constraints on the age of 
stellar populations in SBS 1129+576. 

Emission and absorption hydrogen Balmer lines are 
seen in the spectra of a large part of the galaxy along 
the slit. This allows us to study the age of stellar popu- 
lations, using two methods, based on the time evolution 
of equivalent widths of (1) nebular emission Balmer lines 
and (2) stellar absorption Balmer lines. The advantage 
of these methods is that they are extinction-insensitive. 
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Table 3. Element abundances in 


regions a and b. 


Value 


region a 


region b 


T c (0 ra)(K) 


18930±2540 


20560±2660 


T c (0 ii) (K) 


15350±1930 


lOoOUztZloU 


T C (S m)(K) 


17410±2110 


18760±2210 


l\ c (o llj(cm J 


a i on 
4U±lzU 


10+10 


O+/H+(xl0 5 ) 


1.092±0.351 


1 . Z y O ZC U . O a ^± 


O +2 /H+(xl0 5 ) 


1.215±0.378 


1.691±0.907 


I \ 111 { v in5\ 


z.oU / ztU.Olu 


2.989±1.084 


12 + log(0/H) 


7.363±0.097 


7.476±0.121 


N+/H+(xl0 7 ) 


2.785±0.896 


5.699±2.021 


ICF(N) a 


2.11 


2.30 


log(N/0) 


-1.593±0.241 


-1.414±0.298 


Ne +2 /H+(xl0 5 ) 


0.213±0.070 


0.320±0.165 


ICF(Ne) a 


1.90 


1.77 


log(Ne/0) 


-0.756±0.199 


-0.722±0.245 


He+/H+(A4471) 


0.0686±0.0313 


0.0819±0.0399 


He+/H+(A5876) 


0.0684±0.0094 


0.0603±0.0142 


He+/H+(A6678) 


0.0868±0.0251 


0.0692±0.0447 


He+/H+ 






(weighted mean) 


0.0705±0.0085 


0.0633±0.0129 


He/H 


0.0705±0.0085 


0.0633±0.0129 


Y 


0.2199±0.0271 


0.2019±0.0418 



a ICF is the ionization correction factor. 



This is very important because the only way to derive 
interstellar extinction from optical spectra is to compare 
observed and theoretical decrements of Balmer emission 
lines. However, in the extended low-intensity regions the 
emission lines are weak or absent, making estimates of the 
interstellar extinction uncertain. 

The age of the stellar population in a galaxy can be 
obtained from a third method, comparing the observed 
and theoretical spectral energy distributions (SED), the 
latter computed with various ages and histories of star 
formation for the stellar population. However, the shape 
of the observed continuum is dependent on both age and 
interstellar extinction. If no other information on the stel- 
lar population is available (e.g., the age derived from hy- 
drogen equivalent widths), the extinction coefficient is set 
equal to zero (for a spectrum with no Balmer emission 
lines), and the age can simply be inferred by fitting the 
observed spectrum with theoretical SEDs. This method 
gives a maximum age among possible age estimates. If, 
on the other hand, the ages of stellar populations can be 
estimated by some other methods, then the reddening can 
be derived from the SEDs. 

Finally, to study stellar populations in SBS 1129+576, 
we use broad-band photometry. However, similar to the 
SED fitting method, extinction also affects age determi- 
nation based on photometric data. Additionally, ionized 
gas emission can significantly influence both the observed 
SEDs and broad-band colours. 




5 10 

Distance along the slit (arcsec) 

Fig. 6. Distributions of the Ha equivalent width on a log- 
arithmic scale (thick dashed line), and fluxes of the Ha 
(thick solid line), H/3 (thin solid line), [O in] A5007 (thin 
dashed line) emission lines and the continuum near H/3 
(dotted line) along the major axis of SBS 1129+576. The 
axis origin is set to the position of region a (see Fig. 
The fluxes are in arbitrary units. 



Before dealing with stellar populations, we first con- 
sider how important is the contribution of ionized gas 
emission to the total light. The variations of the equiv- 
alent width of the Ha emission line and the fluxes of 
the strongest emission lines along the major axis of SBS 
1129+576 are shown in Fig. The maximum EW(Ha) 
= 385A and £W(H/3) = 64 A are derived in region a. 
The lower EW(Ha) and EW(R/3) for these regions listed 
in Table [2] are due to different apertures used, of 5" x 2" 
(Table 01 and of 0'.'69x2" (Fig.©. The Ha equivalent 
widths in other regions are much lower. Hence, the contri- 
bution of gaseous emission to broad-band fluxes is small 
in the bright H n regions and negligible in the LSB com- 
ponent. 



5.1. Age from Balmer nebular emission lines 

The largest equivalent widths EW(Ha) are measured in 
regions a and b (Fig.©, implying that the light from these 
regions is dominated by young stellar populations. The 
equivalent widths of Ha, H/3 emission lines in LSB re- 
gions are much smaller, implying a higher contribution to 
the light of the old stellar populations. The fluxes and 
equivalent widths of the Ha, H/3 (and H7 where possible) 
emission lines were measured in the spectra of five LSB 
regions and listed together with errors in Table^J Because 
the H/3 emission line is narrower than the absorption line 
and does not fill the absorption component, its flux was 
measured using the continuum level at the bottom of the 
absorption line. This level has been chosen by visually in- 
terpolating from the absorption line wings to the center 
of the line. 
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The extinction cannot be derived from the Balmcr 
decrement because the observed Ha-to-H/3 emission line 
flux ratios for most of the selected regions are lower than 
theoretical predictions of ionization-boundcd H II region 
models (Tables El 0] and Fig. EJ). Therefore an extinction 
coefficient C(Hfi) of zero is adopted. Two dashed lines 
in Fig. [7| denote the theoretical values F(Ha)/F(Hfi) = 
2.86 for T c = 10 4 K and 2.75 for T = 2 x 10 4 K. Ha 
and Hfi fluxes have been measured in each region except 
for the outermost regions 1 and 5 where only Ha emis- 
sion is present. Only in the two regions, 2 and 4, is the 
Ha/H/3 flux ratio larger than the theoretical recombina- 
tion flux ratio. No correction for the absorption line equiv- 
alent width has been made in these regions. The extinc- 
tion coefficients C(Hfi) for regions 2 and 4 are shown in 
Table 01 Note that CiHfi) for region 4 is very uncertain, 
because of the very low equivalent width EW(Hfi), which 
is comparable to the EW of the Hfi absorption line. 



1 1 

2 


4 


= +4a^-J h 

H 

' 3 ' 

region a 

1 


region b 



-10 10 20 

Distance along the slit (arcsec) 

Fig. 7. The observed Ha/H/3 emission line flux ratios in 
the different regions of SBS 1129+576. Solid horizontal 
bars indicate the apertures wherein the spectra were ex- 
tracted. Large circles show the Ha/H/3 flux ratios for re- 
gions a and 6, computed within different apertures. The 
position of the centers of the two brightest H II regions 
is marked by vertical solid lines. Small circles show the 
same ratios for the fainter regions 2, 3 and 4. The upper 
and lower dashed horizontal lines denote the theoretical 
F (Ha) / F (~H.fi) flux ratios for ionization-bounded (case B) 
Hn regions with T = 10 4 K and T c = 2 x 10 4 K, respec- 
tively. 

The dependence of the Ha emission line equiva- 
lent width on age is shown in Fig. Here we assume 
ionization-bounded Hn region model. The model Ha 
equivalent widths in the case of density-bounded H n re- 
gion are lower than those in Fig. Since the temporal 
evolution of emission line equivalent widths depends on 
the star formation history we consider two limiting cases: 
an instantaneous burst model and models with continuous 
star formation. The instantaneous burst model equivalent 



widths of Ha and Hfi emission lines are calculated using 
the galactic evolution code PEGASE.2 (Fioc & Rocca- 
Volmerange HW7j) . In Fig. Q3 the EW(Ha) for the heavy 
element mass fraction Zq/50 is shown by a solid line and 
that for the heavy element mass fraction Zq/20 by a 
dashed line. PEGASE.2 is based on the Padua stellar evo- 
lutionary models (Bertelli et al. I1994fl and stellar atmo- 
sphere models from Lejeune et al. (|1998fl . An initial mass 
function with a Salpeter slope (a = -2.35), and upper and 
lower mass limits of 120 M Q and 0.1 M Q are adopted for 
all our calculations with the PEGASE.2 code. 

An instantaneous burst model is most appropriate for 
the star-forming regions a and b. The maximum equiva- 
lent widths of Balmer emission lines in region a (EW(Ha) 
= 385 A and EW(Hfi) = 64 A) are consistent with an in- 
stantaneous burst age of 6 Myr for a heavy element mass 
fraction of Z = Z Q /20 and of 9 Myr for Z = Z /5O 
(Fig. [SJ. The maximum value EW(Ha) — 180A for re- 
gion b is consistent with an instantaneous burst age of 9 
Myr for Z = Z Q /20 and 11 Myr for Z = Z /5O. The EWs 
of Balmer emission lines in regions a and b are probably 
slightly underestimated due to the significant (^ 30%; Fig. 
|3Ji) line-of-sight contribution of the emission from the LSB 
component at the respective position. Given the steep de- 
crease of the EW(Ha) as a function of time for 7 Ss log(i 
yr) Ss 7.4 (Fig. [^J this will not significantly reduce the 
instantaneous burst age for regions a and b. 

In continuous star formation models we adopt a con- 
stant star formation rate in the interval between the time 
ii when star formation starts and if when it stops. Time 
is zero now and increases to the past. 

We use the model equivalent widths of hydrogen emis- 
sion lines and SEDs for instantaneous bursts (Fioc & 
Rocca-Volmerange I1997fl to calculate the temporal evo- 
lution of the equivalent widths of hydrogen emission lines 
in the case of continuous star formation with the con- 
stant SFR. The results are given in Fig. |H1 for a heavy 
element mass fraction .Zq/20. The temporal dependence 
of the equivalent width of the Ha emission line is shown 
for continuous star formation starting at time ii, as defined 
by the abscissa value, and stopping at if = 5 Myr (dotted 
line) and if = 8 Myr (dash-dotted line). The equivalent 
width of the Ha emission line in the spectrum of a stellar 
population formed between ii and if is the value of EW 
at time ij (Fig. [SJ. At a fixed EW, the younger are the 
youngest stars, the larger is the time interval ii — if, i.e. 
the older are the oldest stars. 

The models of the continuous star formation apply 
best to the extended stellar component in SBS 1129+576. 
The observed EW(Ha) for regions 1-5 (Table are 
shown in Fig. IHlby triangles on the theoretical curve with 
the star formation stopping at if = 8 Myr ago and corre- 
spond to ages ii in the range 20 - 400 Myr. If instead the 
star formation continues until if = 5 Myr (dotted line), 
then the age i; of the oldest stars, as derived from the 
equivalent widths of the Ha and Hfi emission lines, should 
be in the range 50 Myr - 5 Gyr. However, the latter star 
formation history seems to be less likely if a stellar IMF 
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with upper mass cut-off of M = 100 Mq is assumed. This 
is because massive stars with M = 40 Mq are expected 
to be present in region with if = 5 Myr while the Ha lu- 
minosities of regions 1-5 correspond to fluxes of ionizing 
photons lower than that produced by single massive star 
with M = 40 - 60 Mq. For continuous star formation with 
if = the upper mass cut-off of the stellar IMF must be 
even lower. 




log t, (yr) 

Fig. 8. Time evolution of the equivalent widths of nebular 
emission line Ha for an instantaneous burst with a heavy 
element mass fraction Z = Zq /20 (dashed line) and Z = 
Zq/50 (solid line). Model predictions are also shown for 
the case of constant continuous star formation starting at 
an age defined by the abscissa tj and stopping at if, with 
tf — 5 Myr (dotted line) and if = 8 Myr (dash-dotted 
line), both with Z = Zq/20. Triangles show the positions 
of LSB regions with the observed EW (Ha) superposed on 
the modeled curve for continuous star formation stopping 
at if = 8 Myr. 

The small number of massive stars and hence the 
stochastic nature of the IMF by region can influence 
our age estimates from hydrogen emission line equivalent 
widths. From the observed fluxes of the H/3 emission line 
for regions 1 - 5 we find that 3 - 23 07v or 35 - 233 BOv 
stars are required to ionize the gas (see TableQJ . Here, Lyc 
photon fluxes of 10 49 s _1 and 10 48 s _1 were respectively 
adopted for a single 07v and BOv star (Vacca et al. H99tj[l . 
Then the total stellar mass of the single stellar population 
was estimated to be in the range 10 4 - 10 5 Mq assuming 
an IMF with Salpeter slope and upper and lower stellar 
mass limits of 120 Mq and 2 Mq. The total mass for the 
different regions is in the same range as in the case of SBS 
0940+544 (Guseva et al. l2lM| . Cervino et al. ijMTnjl find 
that in the range of equivalent widths EW(H[3) = 1-10 
A which is typical for regions 1-5 (see Table 0}, the dis- 
persion of age at fixed EW(R(3) is not greater than 5 - 
10% if the total mass of the cluster lies in the range 10 4 - 
10 5 Mq. Hence, we conclude that our age estimate is not 
significantly affected by stochastic effects. 



5.2. Age from the Balmer stellar absorption lines 

We also use a second extinction-insensitive method to es- 
timate stellar population ages from the equivalent widths 
of Balmer absorption lines H7 and H<5. Other higher-order 
hydrogen Balmer absorption lines have not been used for 
age determination due to the relatively noisy spectrum at 
short wavelengths. The ages were derived from the cal- 
ibration of equivalent widths of absorption Balmer lines 
versus age calculated by Gonzalez Delgado et al. (1999) 
for ages 1 Gyr. Their models predict in the case of an 
instantaneous burst of star formation a steady increase of 
the equivalent widths with age for ages ranging from 1 
Myr to 1 Gyr. However, for ages iS 1 Gyr the situation is 
the opposite and the equivalent widths of the absorption 
lines decrease with time (Bica & Alloin 1986). 

The hydrogen absorption lines, created by the under- 
lying stellar populations, are seen in all positions along the 
slit. In two regions, 2 and 3, the H7 absorption line was 
not used for the equivalent width measurements because 
the contamination by the H7 and [O iii]A4363 emission 
lines is too strong (Figs. EH ITT|l . 




7 8 9 10 

log t, (yr) 

Fig. 9. Dependence of the equivalent width of H5 absorp- 
tion line for an instantaneous burst on age for ages > 1 
Gyr from Bica & Alloin (1986) and for ages < 1 Gyr with 
metallicity Z— Zq/20 from Gonzalez Delgado et al. (1999) 
is shown by solid line. Model predictions are also shown 
for the case of constant continuous star formation start- 
ing at an age defined by the abscissa ii and stopping at if, 
with if = 5 Myr (dotted line) and if = 8 Myr (dash-dotted 
line), both with Z — Zq/20. Triangles show the positions 
of the corrected EW(HS) for LSB regions 1-5 superposed 
on the models of continuous star formation with if = 8 
Myr. 

The measured equivalent widths of absorption lines 
were corrected for the contribution of the nebular emis- 
sion in the same way as in Guseva et al. J2Q01). No ex- 
tinction correction has been applied to emission lines su- 
perposed on absorption profiles. Table lists the uncor- 
rected equivalent widths of the H7 and US absorption lines 
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Table 4. Fluxes, equivalent widths of the Ha, H/3 and H7 emission lines and the extinction coefficient C(H/3) in LSB 
regions. 



Region 


Distance 3 


Aperture 13 


F(Hq) c 


EW(Raf 


F(npy 


£W(H/3) d 


F(H 7 ) C 


£W(H 7 ) d 


C(Uf3) 


1 


-7.9 


2.0x6.2 


4.2±0.5 


32.3±2.5 












2 


-5.5 


2.0x5.5 


18.4±0.6 


83.5±1.8 


3.9±0.5 


9.6±0.7 






0.33±0.07 


3 


3.1 


2.0x4.8 


25.4±0.6 


47.3±0.8 


13.5±0.6 


16.2±0.5 


5.6±0.5 


4.2±0.3 




4 


14.1 


2.0x4.8 


4.1±0.3 


23.3±1.9 


1.1±0.3 


2.4±0.7 






0.20±0.06 


5 


17.9 


2.0x4.8 


0.9±0.3 


8.7±2.5 













a distance in arcsec from region a. Negative and positive values correspond to regions located respectively to the southeast and 
northwest from region a. 

b aperture x x y where x is the slit width and y the size along the slit in arcsec. 
c in units 10 -16 erg s -1 cm -2 . 
d in A. 



Table 5. Equivalent widths of H7 and H<5 absorption lines in LSB regions. 



Region 


Distance 81 


Aperture b 


EW{my 


EW(H~f) c 


EW(m) cm c 


£W(H 7 ) cor c 


1 


-7.9 


2.0x6.2 


5.6 ±1.1 


5.3 ±1.3 


5.6 ±1.1 


5.3 ±1.3 


2 


-5.5 


2.0x5.5 


5.7 ±0.6 




6.4 ±0.7 




3 


3.1 


2.0x4.8 


4.7 ±0.4 




7.1 ±0.5 




4 


14.1 


2.0x4.8 


6.9 ±0.5 


4.1 ±0.5 


7.1 ±0.8 


4.6 ±0.5 


5 


17.9 


2.0x4.8 


5.7 ±0.7 


6.5 ±0.9 


5.7 ±0.7 


6.5 ±0.9 



Mistance in arcsec from the brightest H 11 region a. Negative and positive values correspond to regions located respectively to 
southeast and northwest from region a. 

b aperture x x y where x is the slit width and y the size along the slit in arcsec. 
c in A. 



with errors and those corrected for the contribution of the 
nebular emission. The corrected equivalent widths of H7 
and US do not show statistically significant spatial varia- 
tions. Note that the corrections are larger for the brighter 
regions. A careful placement of the continuum is very im- 
portant for deriving accurate EWs. We choose points in 
the spectrum free of nebular and stellar lines for the deter- 
mination of the continuum level. Then the continuum was 
fitted by cubic splines. The uncertainties of the continuum 
level were estimated from several different measurements 
of the equivalent widths of Balmer absorption lines with 
independent continuum fittings. These uncertainties are of 
the same order as the errors in Table [51 obtained from the 
fitting of line profiles with Gaussians and from continuum 
noise. 

The instantaneous burst dependence on the age of the 
equivalent width of the H<5 absorption line (solid line) 
(Bica & Alloin (1986) for ages > 1 Gyr; and Gonzalez 
Delgado et al. (11999) for ages < 1 Gyr at a metallicity Z 
= Z /2Q) is shown in Fig.® 

The temporal evolution of the H7 and US absorption 
line equivalent widths in the case of continuous star forma- 
tion is calculated similarly to that of the Ha and H/3 emis- 
sion line equivalent widths described in Sect. 15.11 More 
specifically, we use the model and empirical equivalent 
widths of hydrogen absorption lines (Gonzalez Delgado et 
al. 119991 Bica & Alloin 11986) and SEDs for instantaneous 
bursts (Fioc & Rocca-Volmerange I1997fl to calculate the 
temporal evolution of the equivalent widths of hydrogen 



absorption lines in the case of continuous star formation 
with constant SFR. The results are given in Fig. |5| The 
temporal dependence of the equivalent width of the US ab- 
sorption line is shown for star formation starting at time 
ti, as defined by the abscissa value, and stopping at tf — 
5 Myr (dotted line) and tf = 8 Myr (dash-dotted line). 



The measured EW(RS) are shown in Fig.[5]on the the- 
oretical curve with the star formation stopping at tf = 8 
Myr ago and are marked by triangles. Their positions are 
consistent with an age ti in the range of ~ 30 - 140 Myr. 
A similar age is obtained from the nebular emission lines. 
If instead the star formation is continuing until t{ = 5 Myr 
(dotted lines), then the age of the oldest stars, as derived 
from the equivalent width of the H<5 absorption line, dif- 
fers significantly from that derived from the Ha emission 
line. Furthermore, the discrepancy between ages derived 
from the equivalent widths of the emission and absorption 
lines is very large for a model with tf = 0. Hence, from 
the comparison of hydrogen emission and absorption line 
equivalent widths we exclude models of star formation in 
which stars are continuously forming with a constant star 
formation rate between and 10 Gyr. Finally, we choose 
the case with tf = 8 Myr to be consistent with the age 
derived from the equivalent widths of the Ha emission 
line. 
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5.3. Age from the spectral energy distribution 

The shape of the spectrum reflects both the properties of 
the stellar population and the reddening. Therefore only 
a combination of the spectral energy distribution method 
with the methods discussed in Sect . 15 . II and 15 . 21 allows us 
to derive simultaneously age and interstellar extinction. 

We use the galactic evolution code PEGASE.2 (Fioc 
& Rocca- Volmerange 2H2Z| to produce a grid of theoreti- 
cal SEDs for an instantaneous burst of star formation and 
ages ranging between and 10 Gyr, and a heavy element 
mass fraction of Z — Z Q /20. Because of the low equiva- 
lent widths of hydrogen emission lines in regions 1-5 the 
gaseous emission is not included in the spectral energy 
distributions. Thus, photometric and spectroscopic data 
give direct information about stellar populations when in- 
terstellar extinction is known. 



5.3.1. Continuous star formation with a young stellar 
population 

To fit the observed SEDs of regions 1-5 with only a young 
stellar population continuously formed over the last 100 
Myr and derive the extinction, we consider star forma- 
tion occurring between t- s — 100 Myr and tf = 8 Myr. 
This model predicts an EW(R6) = -7.2A, EW(JI~f) = - 
5.9k, EW(Rf}) = 3.5k and EW(Ra) = 21.2A, close to 
the values observed in the LSB regions (Tables 0] and EJ) • 
The results of our fitting are shown in Fig. ^| We ad- 
just the extinction coefficient C(H/3) to achieve the best 
agreement between the observed SED, after correction for 
interstellar extinction, and the theoretical SED. The ob- 
served extinction-corrected SEDs are superimposed on the 
synthetic SEDs for regions 1 to 5 in Fig. HOfa — llUb. They 
are labeled by the derived values of C(H/3). The syn- 
thetic SEDs reproduce the observed spectra quite well. 
We therefore conclude that continuous star formation, oc- 
curing during the last 8-100 Myr, is a reasonable model 
of the star formation in the LSB regions. 

We also considered the effect of metallicity on the age 
determination and find that it is small. We fit the observed 
SEDs for regions 1-5 with model SEDs calculated for Z 
= Zq/50 and assuming for each region the same C(H/3) as 
in Fig. ED The best fits are obtained with models charac- 
terised by continuous star formation occuring during the 
last 8-110 Myr, similar to the age range of models with 
the heavy element mass fraction Z = Z Q /20. 

Note, that the extinction, derived from the best fit of 
the observed SEDs, is larger in the southeastern part of 
the galaxy (regions 1 and 2) than in the northwestern one 
(regions 4 and 5). The values for the extinction are similar 
to those in some well studied H n regions, e.g. in the LMC 
(Oey et al. l2000f> . open clusters in our Galaxy (Piatti, Bica 
& ClariaEnnOll and some BCDs (Guseva et al. l307H|) . 




region 5 e) 




4000 ' ' 5000 ' ' 6000 ' ' 7000 
Wavelength (A) 



Fig. 10. Spectra of regions 1 - 5 on which synthetic 
continuum spectral energy distributions are superposed. 
Synthetic SEDs in (a) - (e) are calculated for a stellar 
population formed continuously with a constant star for- 
mation rate between 8 and 100 Myr ago. Extinction coef- 
ficients are derived from the best fits of observed and cal- 
culated SEDs. Each spectrum is labeled by its respective 
C(H/3). The SEDs are superposed on the spectra corrected 
for the extinction. 

5.3.2. Continuous star formation including an old 
stellar population 

We consider next continuous star formation scenarios in 
which an old stellar population is present. For this, we 
adopt C(H/3) = and consider models with constant and 
varying SFRs in the age interval between and 10 Gyr. 
Specifically, for a varying SFR, we consider two periods 
of star formation with constant but different SFRs, oc- 
curring in the age interval < 100 Myr and > 100 Myr. 
To quantify the recent-to-past star formation rate ratio, 
we use the parameter b = SFR(t < 100 Myr)/SFR(i > 
100 Myr) following Guseva et al. (2001). First we consider 
models with constant star formation during the whole 
- 10 Gyr range, i.e. models with 6=1. In Fig. Illb II lb 
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1.0 p — , — . — . — , — , — , — , — , — , — . — . — , — . — , — i — , — ,— 

region 1 a) 



b=7 




region 3 c) 




Wavelength (A) 

Fig. 11. Spectra of regions 1 - 5 on which synthetic 
continuum spectral energy distributions are superposed. 
Synthetic SEDs shown in (a) - (e) and labeled by b = 
SFR(t < 100 Myr)/SFR(i > 100 Myr) = 1 correspond to 
stellar populations forming continuously with a constant 
star formation rate since 10 Gyr ago. Synthetic spectra la- 
beled by b > 1 correspond to stellar populations forming 
continuously between and 10 Gyr with a star formation 
rate enhanced by a factor of b during the last 100 Myr. 
These SEDs are superposed on the observed spectra un- 
corrected for extinction. 

we show such SEDs superimposed on the observed spectra 
of regions 1-5 (labeled b = 1). It is evident that these 
models do not reproduce the observed SEDs. However, by 
increasing the parameter b we can fit the observed SEDs. 
These theoretical SEDs are labeled by b > 1 in Fig. Illb 
- II lb . The predicted equivalent widths of the hydrogen 
emission and absorption lines are also in agreement with 
the observed ones in the case of b > 1. If non- negligible 
extinction is present in regions 1-5, then to fit the ob- 
served SEDs, the parameter b should be further increased. 
In particular, if the extinction coefficient C(H/3) = 0.33 as 
derived for region 2 from the Ha/H/3 flux ratio is adopted, 
then the parameter b should be as high as ~ 100 to fit the 



observations. However, if C(H/3) = 0.33 is assumed for re- 
gion 4, then this model cannot fit the observations, even 
for a b = 1000. This implies, that C(H/3) is likely to be 
small in region 4 (Fig. I10|) . Slightly lower values (by ~ 
10%) of the parameter b are obtained if theoretical SEDs 
with Z = Zq 1 50 are used instead of the models with Z = 
Z Q /20. 

In summary, analysis of the spectroscopic data for 
the LSB regions shows that the stellar population can be 
equally well reproduced by a model of continuous star for- 
mation with constant SFR during the last 100 Myr, or by 
a model in which stars are continuously formed in the pe- 
riod 0-10 Gyr with varying SFRs. In the former case, 
a non-negligible extinction should be taken into account. 
In the latter fit to the observed SED is only pos- 

sible when b > 1 (b = 6 - 50, if C(H/3) = and b > 100, 
if C(H/3) > 0.2 - 0.3), i.e. the star formation rate in the 
LSB component has significantly increased over the last 
100 Myr. 

5.4. Age from the colour distribution 

Photometric data allow us to check the results obtained 
in Sect. 5.1 through 5.3. For this purpose we extract from 
the V and I images the same areas as those covered by the 
long slit spectroscopic observations and compare the (V — 
I) colours with predictions from our population synthesis 
modeling. 

The results of this comparison are shown in Fig. 
With solid lines we show in panels (a) and (b) the observed 
V and I surface brightness distributions in mag arcsec -2 
along the slit while the corresponding (V — I) colour dis- 
tribution is shown in panel (c). The predicted colours, ob- 
tained from convolved theoretical spectral energy distribu- 
tions with the appropriate filter bandpasses are shown in 
Fig. 112b with different symbols. The transmission curves 
for Johnson V and Cousins I bands are taken from Bessell 
(1990) and the zeropoints are from Bessell et al. |1998). 

Since the contribution of the gaseous emission is very 
small in regions 1-5, we do not take it into account and 
consider the colours of a stellar population formed contin- 
uously with a constant star formation rate between 8 Myr 
and 100 Myr ago and reddened by interstellar extinction. 
In Fig. 112b we show by small filled circles the colours of 
five regions with extinction coefficients derived from the 
best theoretical fit to the observed extinction-corrected 
spectra (see spectra in Fig. Milk -HTIfeV Open circles show 
the predicted colours of a stellar population continuously 
formed between and 10 Gyr ago, assuming an enhanced 
star formation during the last 100 Myr, as defined by the 
parameter b (see upper spectra in Fig. II lh -lllfc). In this 
case C(H/3) = is adopted. 

The contribution of the ionized gas emission to the 
integrated light of regions a and b is small but not neg- 
ligible. We have, therefore, included gaseous emission in 
our calculations. We have constructed stellar population 
SEDs for a heavy element mass fraction Z = Z Q /20 and 
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Fig. 12. (a) and (b) Surface brightness distributions along 
the slit in V and I. The origin is at the location of region 
a. (c) (V ~ I) colour distribution along the slit. The two 
bright H n regions are marked as "region a" and "region 
b" . The LSB regions are marked as "1" to "5" . Small filled 
circles show the predicted colours of a stellar population 
formed continuously with a constant star formation rate 
between 8 Myr and 100 Myr ago and reddened with an 
extinction coefficient as derived for 5 regions (see Fig. ll(J|) . 
Open circles indicate the predicted colours of stellar pop- 
ulations formed continuously between and 10 Gyr with 
the enhancement parameters b as derived for these regions 
(see Fig. lll(l . Large filled circles show the colours for a stel- 
lar population formed in an instantaneous burst with ages 
t = 8 and 10 Myr for regions a and b, respectively. These 
ages are derived from the observed equivalent widths of 
the Ha and H/3 emission lines. 



instantaneous burst ages of 8 and 10 Myr for regions a 
and b respectively as derived from EW(Ra) (Fig. [HJ and 
EW(Hf3). The calculation of the gaseous SED is made in 
the same way as in Guseva et al. (2001) and then added to 
the SED of the stellar component. The derived colours are 
shown in Fig. 112b by large filled circles. The colour varies 
sharply near the peak of intensity in regions a and b. The 
agreement between the observed and calculated colours in 
the wide apertures is fair, whereas it is better in the small 
aperture of region a. 



The agreement between the (V — I) colours obtained 
from the photometric data and those derived from the 
spectral energy distributions in both scenarios is quite 
good, given the large extraction apertures for regions 
along the slit and a possible small positional shift between 
photometric and spectroscopic data. The systematically 
redder colours observed at the relative position of regions 
4 and 5 are likely caused by the red background galaxies 
such as bg (Fig.QJ, intersecting the slit in the northwestern 
part of the galaxy. This emphasizes the need for identifi- 
cation and rejection of background or foreground sources 
when colours are used to derive ages. Figure 1 shows a 
sizable number of red background galaxies all over the 
field of SBS 1129+576 spanning a range of about 10 mag 
with respect to their I magnitudes. Therefore, the failure 
of any continuous star formation model to reproduce the 
observed red colours in regions 4 and 5 is not surprising. 
Note that the average (V— I) colour in the outer part of 
SBS 1129+576 does not exceed 0.6 mag, i.e. it is bluer by 
^ 0.2 mag than the observed colour in regions 4 and 5. 

6. Conclusions 

We present a detailed spectroscopic and photometric 
study of the very metal-deficient dwarf irregular galaxy 
SBS 1129+576, a likely young galaxy candidate. Broad- 
band V and / imaging and spectra in the optical range 
have been obtained with the Kitt Peak 2.1m and 4m tele- 
scopes, respectively. The main conclusions of this study 
can be summarized as follows: 

1. SBS 1129+576 is a very low-metallicity nearby (D 
= 26.3 Mpc) dwarf galaxy with a chain of H n re- 
gions along its elongated stellar low-surface brightness 
(LSB) component. The average (V — I) colour of the 
outer part of the galaxy with a surface brightness \iy 
in the range 23 - 26 mag arcsec -2 is relatively blue 
~0. 56+0. 03 mag, as compared to (V — I) ~0. 9-1.0 in 
the majority of dwarf irregular and blue compact dwarf 
(BCD) galaxies. The scale length a obtained from the 
surface brightness profiles in the V and / is ~ 430 pc. 

2. The oxygen abundance is found to be respectively 12 
+ log(0/H) = 7.36 ± 0.10 and 7.48 ± 0.12 in the two 
brightest H II regions (regions a and b). Because of the 
low intensity of the [O iii]A4363 in region a, shock en- 
hancement of this line may be important, which may 
result in a slight underestimate of the oxygen abun- 
dance. The neon-to-oxygen abundance ratio log Ne/O 
= -0.76 in region a is in good agreement with the 
mean ratio derived from the previous studies of low- 
metallicity galaxies (e.g., Izotov & Thuan 1999). The 
nitrogen-to-oxygen abundance ratio log N/O = -1.60 
lies in the narrow range of the N/O ratios obtained by 
Thuan et al. ifTMS) and Izotov & Thuan fTUUty for the 
most metal-deficient BCDs. 

3. The He mass fractions Y = 0.220 ± 0.027 and 0.202 
± 0.042 derived, respectively, in regions a and b are 
significantly lower than the value of the primordial He 
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mass fraction Y p = 0.244 - 0.245 derived previously by 
Izotov & Thuan (1998) and Izotov et al. (JIHUJ). This 
difference is likely due to significant underlying He I 
stellar absorption in SBS 1129+576. Hence, despite its 
low metallicity, this galaxy is not a good candidate for 
primordial helium determination. 

4. Hydrogen Ha and H/3 Balmer lines are seen in emis- 
sion in the LSB component, while higher order hydro- 
gen lines are in absorption. Two extinction-insensitive 
methods, based on the temporal evolution of the Ha 
and H/3 emission line and the H7 and US absorption 
line equivalent widths are used for age determination. 
A third method (not extinction-independent) is based 
on the age dependence of the spectral energy distribu- 
tion. Several star formation histories have been consid- 
ered. We find that models of continuous star formation 
with a constant star formation rate starting 10 Gyr ago 
are excluded. However, models starting 10 Gyr ago and 
continuing to the present with a varying star formation 
rate are able to account for the observed properties of 
the LSB regions. Models with star formation rates en- 
hanced by 6 - 50 times during the last 100 Myr can 
reproduce the observed equivalent widths of the emis- 
sion and absorption hydrogen lines and SEDs if zero 
extinction is assumed. If some extinction is present in 
the LSB component then the star formation rate dur- 
ing the last 100 Myr should be enhanced by a factor 
of more than 100 times, or alternatively the observed 
spectroscopic and photometric characteristics of the 
LSB component can be reproduced by models in which 
stars were continuously formed during the last 100 Myr 
only. 

5. The observed (V — I) colours in the LSB component 
of the galaxy are consistent with colours of synthetic 
SEDs for all above mentioned histories of star for- 
mation. Hence, we conclude, that there is no com- 
pelling evidence for either a young or an old age of 
SBS 1129+576. 
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